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Abstract 

We study information geometry of the thermodynamics of first and second 
order phase transitions, and beyond criticahty, in magnetic and hquid sys- 
tems. We estabhsh a universal microscopic characterization of such phase 
transitions via the equahty of correlation lengths ^ in coexisting phases, 
where ^ is related to the scalar curvature of the equilibrium thermodynamic 
state space. The 1-D Ising model, and the mean-field Curie- Weiss model 
are discussed, and we show that information geometry correctly describes 
the phase behavior for the latter. The Widom lines for these systems are 
also established. We further study a simple model for the thermodynamics 
of liquid-liquid phase co-existence, and show that our method provides a 
simple and direct way to obtain its phase behavior and the locations of the 
Widom lines. Our analysis points towards multiple Widom lines in liquid 
systems. 
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1 Introduction 



The physics of phase transitions has been a fascinating area of research for 
more than a century, starting from the celebrated van der Waals equation for 
hquid-gas systems. Apphcation of information geometric methods to such stud- 
ies are relatively recent, but have yielded several useful insights especially in 
the context of liquid-gas co-existence |2J. This method involves the underlying 
Riemannian geometry of the equilibrium thermodynamic state space, and was 
mainly initiated through the work of Weinhold [3] and Ruppeiner [3] |1| The key 
idea here is to utilise the positivity condition on the Hessian of the entropy for 
a thermodynamic system in equilibrium, so as to define a Riemannian metric on 
the thermodynamic state space. In a striking conjecture, Ruppeiner proposed, 
via the theory of Gaussian fluctuations, that the scalar curvature R (or, more 
appropriately, \R\), arising out of such a metric is related to the correlation length 
^ of the system, an idea that has since been tested in a variety of models. Indeed, 
it has been established by now that the scalar curvature of the thermodynamic 
state space diverges at critical points, for a wide variety of systems that exhibit 
second order phase transitionso 

Most of the analysis done in the past dealt with the geometry of thermody- 
namics at or near criticality. However, it has recently been established [6] that 
methods of information geometry can be used to study discontinuous first order 
phase transitions, via the equality of the correlation lengths of the coexisting 
phases, in liquid-gas systems. In [6], such equality has, in fact, been verified with 
experimental data from NIST [7]. Apart from providing an alternative route 
to study phase transitions that bypasses several long-standing problems in stan- 
dard thermodynamics, this approach has the potential of analytically predicting 
the "Widom hne" [8], [9], a conjectured continuation of phase co-existence that 
extends beyond criticality, and distinguishes between supercritical "phases" of 
a system. The presence of the Widom line, which is usually defined as the lo- 
cus of maxima of the correlation length has recently been experimentally estab- 
lished [To], and certainly puts information theoretic studies of phase transitions 
in perspective. Note that very close to criticality, all response coefficients scale as 
powers of the correlation length, and hence the loci of maxima of these serve as 
equivalent definitions of the Widom line [8]. Slightly away from criticality, how- 
ever, this is not the case and our method provides a way of locating the Widom 
line as per its actual definition, even away from criticality. 

Whereas methods of information geometry in liquid-gas systems have been 
firmly established by now, magnetic and liquid-liquid systems have been much 
less studied. Although it was known that the thermodynamic scalar curvature 
can be calculated for some magnetic systems, and these show the standard di- 
vergences associated with a second order critical point, it is not known how 
Riemannian geometry captures first order transitions in these systems. Even 
less is known about the information geometry of liquid systems and its role in 

^Information geometry in the context of thermodynamics is also cahed "thermodynamic 
geometry." 

^Information geometry of quantum phase transitions has also been an area of extensive 
research, see e.g. |5j. 
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predicting phase behavior therein. Here we address these issues and show that 
geometrical techniques can be effectively used to study phase behavior and the 
Widom line in these cases as well. Our results here, in conjunction with those 
for liquid-gas systems analyzed in [6] firmly establish the usefulness and applica- 
bility of information geometry in mean field thermodynamic systems, on which 
experimental predictions can be based. We note here however, that although the 
equality of correlation lengths in co-exisiting phases in liquid-gas systems can 
be understood [6] following the work of Widom |Tl], such a picture has not be 
rigorously developed for magnetic and liquid systems. We will proceed by as- 
suming this, and show that it exactly reproduces known features in mean-field 
magnetic models. With this understanding, we then analyse liquid-liquid sys- 
tems in a simplified mean-field model. The usefulness of our method is that it 
provides an analytic tool for predicting phase transitions and the Widom line in 
difficult situations where one otherwise needs to resort to molecular dynamics 
simulations [9]. 

This paper is organised as follows. In the first section, we briefly review the 
information geometry of equilibrium thermodynamics, and then proceed to study 
the exact solution of the 1-D Ising model in this context. This section is mostly 
devoted to establishing the notations and conventions that we follow in the rest of 
the paper, but we present an interesting observation regarding the Widom line in 
exactly solvable systems. Next, we move on to study phase behavior of magnetic 
systems in the mean-field Curie- Weiss model, and finally analyse liquid systems 
in section 3. For both these cases, the prediction of phase co-existence and the 
Widom line is elaborated using geometric methods. We end the paper with our 
conclusions and some future prospects. 

2 Information Geometry and The 1-D Ising Model 

We begin with a very brief overview of information geometry. It will be enough for 
us to consider this in the context of equilibrium thermodynamics, a formulation 
due to Ruppeiner |4j. More details can be found in the excellent review of Brody 
and Hook [2J. The main idea is to consider the line element due to a positive 
definite Riemannian metric which is defined by the Hessian of the entropy per 
unit volume s, 



where and denote the internal energy and the particle number per unit vol- 
ume (for magnetic systems, one has to appropriately use thermodynamic quanti- 
ties per unit spin). As shown by Ruppeiner [1], this introduces the concept of a 
distance in the space of equilibrium thermodynamic states, i.e, a large distance 
between two such equilibrium states is interpreted as a small probability that 
these are related by a thermal fluctuation. 

One can consider various other forms of the metric related to the one given 
in eq. ([1]) by Legendre transforms. Indeed, a choice of the thermodynamic poten- 
tial in appropriate coordinates might render the metric diagonal, and simplify 
algebraic calculations. For example, for single component fluids and magnetic 




(1) 
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systems, a particularly simple diagonal form of the metric can be used 



kBT\dTj^ kBT\dpJ^ 

where fi = (^^^ y f being the Helmholtz free energy per unit volume (per unit 

spin for magnetic systems), and p is the inverse of the volume. For magnetic 
systems, the magnetization per unit volume plays the role of p. 

For such diagonal metrics, the scalar curvature of the manifold takes a well 
known simple form. 
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d f 1 dgpp\ ^ d f 1 dgTT\ 



dT\^ dT J dp \^ dp J 



(3) 



where qtt and Qpp are the coefficients multiplying dT"^ and dp^ respectively, in 
eq.([2]), and g = gxT-Qpp- In passing, we note that for positivity of the line element, 
Qtt and Qpp should be positive definite in the regions of interest. Non-diagonal 
forms of the metric (as in eq.([T])) can also be used, and give equivalent results for 
R. 

By augmenting this geometric picture with Gaussian fiuctuation theory, Rup- 
peiner conjectured that such a scalar curvature can, in fact, be related to the 
correlation length of the system \R\ ~ where d is the system dimension. 
We also mention that the geometry described here becomes trivial in the presence 
of a single fiuctuating variable, and meaningful results can only be obtained for 
two or more fiuctuating thermodynamic quantities. 

In the context of the 1-D Ising model, information geometry has been studied 
by Ruppeiner ^12j and later by Janyszek and Mrugala, [13]. We briefiy recall 
their results. We start with the thermodynamic potential per unit spin in the 
ferromagnetic case (with unit coupling constant), which is given by 



$ = -^In 



e^cosh (a) + (e'^sinh^ (a) + e"'^) '^'1 (4) 



where H and T denote the applied magnetic field and the temperature respec- 
tively, and a = H/T, P = 1/T (the Boltzmann's constant has been set to be 
unity). From eq.dl]), the authors of [13] showed that in the (a,/3) coordinates, 
the scalar curvature of the thermodynamic state space is given by the simple 
analytical expression 

':o^^(o:) (5) 

(smli^(a)+e-«)'" 

This also agrees with the result of [12]. In the zero field limit, the curvature 
diverges at T = indicating the sole critical point of this theory. Also, it can be 
checked that 

^'^-°= |ln(tanh(/3))r ^''^'^-° 
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Figure 1: "Widom lines" for the 1-D Ising model on the H — T plane. The red 
curve is the locus of maxima of the scalar curvature R. The blue and the green 
curves are the maxima of the specific heats Cm and Ch respectively. 



In order to understand the Widom line, we study information geometry of the 
1-D Ising model away from criticality for non-zero T and H. We compute the 
maxima of the correlation length (equivalently, |-R|) in the supercritical region, 
and contrast it with the maxima of other response functions. We first record the 
standard expression for the magnetization per spin 

sinh (a) , , 

" (sinh^ (a) + e-4/3) 

The specific heat at constant magnetization is then obtained by differentiating 
the entropy obtained from eq.Q, while keeping m fixed, 

4p2a 

r 



T2 (e^/^sinh^ (a) + l)'/' (^e^/^cosh (a) + (e^/^sinh' (a) + l)'^') cosh 



a 



with a lengthy expression for the specific heat Ch at constant applied field H 
which is not reproduced here. These of course reduce to the standard expression 
Cm,H\H=o = /3^sech^/3 in the zero field case. Now, we calculate the maxima of 
\R\ as a function of temperature. From eq.(l5]), it can be seen that such maxima 
occur for a given value of H when the the temperature T satisfies the equation 

and a similar analysis can be carried out to determine the maxima of the specific 
heats Cm and Ch- El Numerical solutions of these have been plotted in fig.([T]), 



We will loosely refer to the locus of maxima of the correlation length and the other response 
coefficients as the "Widom lines." 
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where the red, green and blue curves denote the locus of maxima of \R\, Ch, and 
Cm respectively. We observe that near the critical point, i.e for very small values 
of the applied magnetic field H, the maxima of the specific heat do not asymptote 
to T = 0, but reaches a limiting temperature T ~ 0.8. The maxima of \R\ (i.e 
the correlation length ^) for such small values of H also does not asymptote to 
zero, but goes to a limit T ~ 0.18 for H = 10~^. It is interesting to note that the 
locus of maxima of the correlation length does not begin from the critical point 
for the 1-D Ising model. As we will see, this feature will not be present in any of 
the mean-field models that we will analyse. 

Before concluding this section, let us summarize the main results here. We 
have considered the 1-D Ising model from the point of view of information ge- 
ometry. Analysing of the same away from the critical point, we have located the 
locus of maxima of the correlation length ^ and the specific heats Cm and Ch, 
in the T — H plane. Our analysis points to the fact that in this simple example, 
the Widom lines associated with the locus of maxima of R and also the specific 
heats do not seem to converge to the critical point (T= 0) in the limit H ^ 0. 



3 Information Geometry and The Curie- Weiss 
Mean Field Model 

In this section, we study the information geometry of the classical mean-field 
Curie- Weiss (CW) ferromagnetic model. This model has been studied extensively 
in the past, and details can be found in standard textbooks [H], [15]. We will 
adhere to the notations of [T^. In the context of information geometry, the 
CW model was first studied by Janyszek and Mrugala in [13\. These authors 
established the divergence of the scalar curvature close to the critical point. We 
seek to understand this model away from criticality, and study first order phase 
transitions, via the geometry of the thermodynamic state space. 

As is well known, in the thermodynamic limit, the free energy of the CW 
model is given by 

= -T(ln2 + max/(m)) (10) 

where the Boltzmann constant has been set to unity, and m is the magnetization 
per spin, which solves the equation 

m = tanh —m + — (11) 

Here, T is the temperature, Tc its critical value, and H is the applied magnetic 
field. In order to write this in terms of the temperature and the thermodynamic 
extensive variable, m, we perform a Legendre transform and write the free energy 
as [13] 

1 T 

/ = -Tln2 - ^T^m"^ + -In (l - m^) + Tmtanh~^m (12) 

with {§^)rp = H. This can be used to calculate the scalar curvature of the 
equilibrium thermodynamic state space, but note that the entropy per spin is 
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given by eq. (fT2|) as 



s = ln2 — -In (l — m^) — mtanh (13) 

As can be immediately seen from eq.f llSp and eq.(IT]), information geometry be- 
comes trivial in this context, since the entropy is a function of a single variable, 
and is not amenable to a Riemannian geometric analysis. In order to remedy the 
situation, the authors of [T3] proposed to modify the Hamiltonian, which in its 
original form for N spins, is given by 

n = -'^J2'^'^-^J2'i^ hj = i---N (14) 

i<j i 

by an additional term that corresponds to the mechanical energy of the lattice. In 
particular, the modification seeks to retain a non zero specific heat for the model, 
which is otherwise zero as can be seen from eq.( !T3|) . Doing this modification ad 
hoc essentially amounts to assuming that 

-^(S) =Cl{T) (15) 

\ / m 

where Cl, the "lattice specific heat" (not to be confused with the specific heat 
at constant m) is an apriori unknown function of the temperature. As we show 
below, this modification nevertheless captures the correct information theoretic 
properties of the system. 

The line element is now given by [13] 

where we have used the free energy representation, with coordinates (T, m) |4j . 
A simple calculation then shows that the scalar curvature for this model is given 
hjR = A.B, where 

^ 1 (1 - m^) 

2C2(T-r,(l-m2))' 
B = T,m\TC'^ + CL)-T,{CL + TC'j} + T{2CL + TC'L) (17) 

Here, the primes denote derivatives of Cl with respect to the temperature T. In 
particular, we see that the curvature diverges as T —)■ Tc (i.e m — )■ 0), as derived 
in a slightly different way in [13]. We now present a general analysis of i?, valid 
for non-zero T and H, away from criticality. The form of the curvature scalar in 
eq.(|T71) is useful when we assume a specific form of Cl and a value of Tc. For 

''Note that the second term in eq.(fT6|) is also related to the inverse of the magnetic suscep- 
tibihty, i.e [B)=^ = ^ 
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Figure 2: Isothermal R vs m for the 
Curie-Weiss model, below criticality. 
The red, magenta and green curves cor- 
respond to T = 0.4, 0.8 and 0.9 respec- 
tively. The critical temperature has been 
chosen to be = 1. 



Figure 3: Isothermal R vs m for the 
Curie-Weiss model, beyond criticality. 
The green line shows a diverging curva- 
ture at the origin for T = Tc = 1. The 
magenta and red curves correspond to 
T = 1.1 and 1.2 respectively. 



simplicity, we assume Tc = 1 and a power law, Cl (T) = 1 + T + T^. This is only 
for illustration, and as can be checked, any other value of the critical temperature 
or any other regular functional form of will not alter our discussion below. 
With these assumptions, the curvature simplifies to 



^ 1 (1 - m2) (H- 2T + 3T2) + 4T3 - 1] 

R = 7i o ( lo ) 

2 ^i^m^ ^Tf{l + T + T^f ^ ' 



We also calculate the specific heat at constant applied field if, using 

H = Ttanh~^m — TcTti (19) 
and obtain from eq.f lT^ . for = 1, 

r r(^'\ T(tanh-^m)^(m^-l) 

We observe that the curvature scalar blows up wherever Ch does, as expected, 
for T < Tc. (However, Ch does not diverge at the critical point). Further, the 
denominator of R behaves as the product of the square of and the square of 
the denominator of Ch- We keep this in mind, as this seems to be a universal 
feature in mean-field theories. 

Now, we study the behavior of i? as a function of m. This is shown in fig. ([2]) 
for temperatures below criticality and in fig. ([3]) for temperatures at and beyond 
criticality. For T < T^, we find that R diverges symmetrically on the m-axis for 
a given value of temperature. This value of m can, in fact, be identified with 
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Figure 4: R ys H for isotherms the 
Curie-Weiss model. The green curves 
correspond to \R\ in the two physical re- 
gions at T = 0.6, and the blue curves 
are for T = 0.8. The red curve corre- 
sponds to the isotherm T = 1.35. Cross- 
ing of the physical branches of indi- 
cate phase transition, always ai H = 0. 



Figure 5: Loci of maxima for the CW 
model. The green line along the H- 
axis denotes phase co-existence and ter- 
minates at Tc = 1. The red line, along 
H = are the maxima of |-R| (and Xt), 
the Widom line for various values of T. 
The two blue lines are the symmetric 
maxima of Ch as a function of T. 



the turning point of the magnetic isotherms in the H — m plane. We have not 
shown these isotherms here, but we note that for values of m that lies between 
the two divergences, the specific heat Ch is negative, i.e, the system becoms 
unstable. The two divergence on the m-axis get closer as we approach Tc = 1, at 
which point they merge into a single divergence at m = 0. Beyond criticality, the 
situation is depicted in fig.(|3]). Here we see that isothermal R has a maximum 
at m = 0, with its value at the maximum decreasing as we move away from 
criticality. 

It is more useful to consider the behavior of i? as a function of the applied 
magnetic field, if, with m being treated as a parameter. We can guess the result 
by exploiting the symmetry of the situation. From the expression for R of eq. ffTSj) . 
we see that R (m) = R (— m). For two physical branches (where Ch is positive), 
the i?'s should thus "cross" where m = 0. From eq. flT^ . this occurs when H = 0. 

We show this in fig.(|l]), for the isotherms T = 0.6 (green), T = 0.8 (blue) 
and T = 1.35 (red), which we now explain this in details. Take, for example, the 
isotherm T = 0.8 of fig. ([2]) (the magenta curve). For this isotherm, the two blue 
curves of fig.(jlD are plots of R in the the regions where Ch is positive, with the 
curve that asymptotes to infinity for positive values of H are for negative values 
of m and the other one is for corresponding positive values of m. Specifically, 
from fig.(l2]), the T = 0.8 i?-isotherm diverges at m = ±0.4472 (as can be seen 
from eq.f jTSj) . The blue curve in fig.dl]) that diverges at a positive value of if is a 
parametric plot for m < —0.4472 and the one that diverges at a negative value 
of H is for m > 0.4472. From eq.( fT9l) . the divergence occurs at if = ±0.0622. 
This is exactly as shown in fig.(j2]). Beyond Tc, there is no i?-crossing. 
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We thus see that below the critical temperature, R has two physical branches, 
which cross at H = 0. Indeed, as alluded to in the introduction, this implies 
equality of correlation lengths at if = which we interpret as the phase tran- 
sition, with the residual magnetization being the values of m where R diverges 
for this temperature. Of course beyond Tc, there is no residual magnetization, as 
can be seen from the behavior of R in fig. ([3]). This is the information geomet- 
ric description of discontinuous transitions in magnetic systems, and agrees with 
standard results. 

As we move closer to Tc, the value of R at which the physical branches cross, 
increases. This can be seen by comparing the green and blue curves in fig.dl]), 
which denote physical branches of i? at T = 0.6 and T = 0.8 respectively. At Tc = 
1, the crossing point is pushed to infinity. Beyond Tc, R shows a maxima, which 
is always at H = 0, which can again be justified using symmetry arguments. The 
location of the Widom line is hence along the if-axis, evan away from criticality. 
This is depicted in fig.(j5]), where the green line {H=0), is the first order line that 
culminates at Tc = 1, and the locus of maxima of R is its continuation, along the 
red line. 

To contrast the behavior of the maxima of R with other response coefficients, 
we have also plotted, in fig.©, the locus of maxima of the specific heat Ch- We 
find that unlike R (and the susceptibility xt, whose locus of maxima is the same 
as that of R), Ch shows two symmetric maxima for a given temperature on the 
if-axis. All maxima converge to the critical point for very small values of H. 

To conclude this section, we summarize the main results. Here, we have 
provided a novel characterisation of first order processes in a simple mean-field 
magnetic model, the CW model of ferromagnetism, via information geometry. 
We have seen that equality of the correlation length ^ of co-existing phases accu- 
rately predicts known behavior of this model, near or away from criticality. The 
interpretation of the scalar curvature of the thermodynamic state space R further 
allows us to calculate the Widom line as the locus of maxima of the correlation 
length C, (via R ~ C,"^) which is shown to lie on the H-axis, as an expected contin- 
uation of the phase co-existence line, in the T — H plane. Such loci of maxima 
for the specific heat Ch however do not lie on the if-axis. We further note that 
our information geometric study of first order processes in the CW model has 
essentially similar features as mean- field liquid-gas systems [6], as expected, and 
indicate the wide applicability of this method in studying phase transitions. As 
an illustration of this, we now proceed to study liquid systems. 

4 Information Geometry and Liquid Systems 

We now apply our information geometric method to a simplified toy example in 
the thermodynamics of liquid-liquid phase transitions. Although the literature 
on the topic is vast (see, for example, [16] for a recent analysis in silica) analytic 
insight in such systems is often difficult to provide in practice. We choose a simple 
model that nevertheless captures the essential physics. The model that we study 
has internal energy 

U = A{V) + B {V) T^/s + kT (21) 
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Figure 6: Isotherms in the V — P plane Figure 7: Meta-stabihty region in the 
for various values of temperature. T — V plane. 

where k is related to Boltzmann's constant. This form of the energy is motivated 
from the fact that the potential energy U ~ for some simple liquids [T7] . 
|18] . Here, A {V) and B {V) are functions of the volume V, which can be taken 
to be polynomial fits, determined from simulation data. The exact nature of 
these functions for liquid silica has been dealt with extensively in [18j (where 
a similar form of the internal energy was used), and they were fitted to fourth 
order polynomials in the volume. For analytical tractability, it is enough for us 
to choose a simpler situation, and we assume 

A(y) = Ao + AiV + A2V^ 

B{V) = Bo + B^V (22) 

where Ai, i = 0---2 and Bj, j = 0,1 are constants. Inclusion of a quadritic 
term in B [V) does not allow for analytic handling of the model, and we have 
dropped this for the time being. Indeed, this simple choice of the functions A and 
B will serve to illustrate the main features for any physical system with a T^/^ 
dependence of the internal energy, and it can be checked that the inclusion of a 
quadritic term in B (V) will not alter the qualitative aspects of our discussion. 
We will also assume throughout that the units are properly chosen, as in [T8] . 
Given the internal energy U, one can compute the entropy of the system by text- 
book methods ^J. In particular, starting from a reference entropy S (Tq, Vq) at a 
reference temperature Tq, and a reference volume Vq, the entropy at arbitrary val- 
ues of the temperature and volume can be calculated by adding its change along 
and isotherm and along an isochore. This requires the expression for the pressure 
for the reference temperature, which can be again obtained as a polynomial fit 
in the density [18]. In our example, we choose 

P\t=To = Co + y + ^ (23) 

Where the Cs are coefficients that can be calculated from molecular dynamics 
simulations. The calculation for the entropy has been pedagogically explained 
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in [18], and we obtain 

S = ^^^[T'/^(2V\nT + 5VBiV) + 

+ 2VA (V) + 2{V + CoV"^ + CiVlnV - C2) ) - 3VB (V)] (24) 

where we have set the constant k in eq.f l2T|) . and the reference temperature Tq to 
unity, without loss of generahty. The Helmholtz free energy for this model can 
be calculated from F = U — TS, and taking its derivative yields the pressure 

P = TPo - (T^/^ - T) + (T - 1) (Ai + 2^2 V^) (25) 
The curve of metastability can be obtained by setting 

^ = 2^2 (T - 1) - {C,V + 20 = (26) 

The critical volume is calculated by solving for T from the above equation and 
then using (^^) = 0, where Tm is the solution for the temperature obtained 
from eq. (l26l) . These finally yield the critical temperature and volume : 

T-^^^^ V--3^ (27) 

The constant volume specific heat for this model can be calculated to be 

5r^/5 + 3(Po + PiV) 

= (28) 

And the isothermal compressibility is given by 

1 (dV\ 

^ ~V \dP)^ ^ T {CiV + 2C2) - 2A2V^ (T-l) ^^^^ 

The expression for Cp, the specific heat at constant pressure is easy to obtain 
by differentiating eq.(l2^ while keeping the pressure in eq.( l25i) fixed, but gives a 
lengthy expressions, which we shall not present fully here, but we note here that 
it has the same denominator as eq. (l29l) . as expected. 

Our main interest here is in the scalar curvature of the equilibrium thermo- 
dynamic state space of this system. This can be calculated with the help of 
either eq.([l]) or eq.(l2]), and we obtain an expression for R in terms of T and V. 
Using eq.(l2]) for example, we obtain a diagonal form of the metric, given by the 
components 

9TT - 9pp - (30) 

where p = 1/V. As required these are certainly positive in the domain of 
interest, as is the quantity ^JQttQpp- The expression for R is too lengthy to 
reproduce here, but for the moment, let us note that its denominator is 

Pden = [(ST^/s + 3 (Po + B^r)) (T (CiV + 2C2) - 2A2V^ {T - 1))]' (31) 
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Figure 8: i? as a function of volume for 
various isotherms. The red, green and 
blue curves correspond to T = 1.6, 1.8 
and 2.1 respectively. Beyond criticality, 
\R\ shows two local maxima. 
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Figure 9: Maxima of |-R| (black and 
brown curves), Cp (red) and Kt (blue) 
in the V — T plane, beginning from 
the critical point in the spinodal (green) 
curve in the V — T plane. 



implying that the curvature scalar diverges along the divergence of Kt or Cp. 
Also note that the denominator of R is equal to the product of the square of the 
numerator of Cy and the square of the denominator of Cp (or Kt)- This is a 
feature that we noticed in the Curie- Weiss model of the previous section, and 
seems to be an universal property for mean-field systems. 

In order to make our results more tractable in general, we now make a choice 
of constants, and set A2 = —1, Bq = 5 x 10^, Cq = 10, Ci = —3, with all the 
other constants set to unity. One can check that this ensures the positivity of 
physical quantities like temperature, volume, isothermal compressibility etc. in 
some domain, with our analysis being valid in this domain. With this choice of 
parameters, the metastability condition of eq.(l26ll yields, as a solution for the 
temperature, 

= 2^3 _ 3V + 2 ^^^^ 

with the system being stable outside the region defined by the above equation, 
plotted in the T—V plane. The critical temperature for our model is, from eq. (l27|) . 
Tc = 2, with Vc = 1 and Pc = 16.21 (in appropriate units). In fig.(|6]), we have 
shown some isotherms on the P — V plane, along with the region of metastability 
for our model. The latter region is plotted in the T — V plane in fig. ([7]). We 
now present our results on the scalar curvature graphically. We find that below 
Tc, isothermal R, as a function of the volume, diverges at the boundaries of the 
metastability curve of fig. ([7]). This is shown in fig. ([8]). Consider, for example, 
the isotherm at T = 1.6, plotted in red in fig. ([8]). For this temperature, Cp 
diverges for V = 0.7899 and V = 1.4844. These are also the values of V for 
which R diverges, and as can be seen from fig.(I7]) (or, equivalently, from eq.( l26l) 
after putting in the constants), these are the values of V between which the 
system becomes unstable. As we look at isotherms with higher temperature, for 
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Figure 10: Isothermal R-crossing as a 
function of P. Two curves of the same 
color denote \R\ in the physical regions 
for a given isotherm. The red curves 
correspond to T = 1.6, the green to 
T = 1.65 and the blue to T = 1.7. The 
physical regions of the P — V isotherms 
have been used for the plots. 
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Figure 11: Widom lines for the model 
liquid-liquid system. The green and 
black curves correspond to the two max- 
ima of |-R|. The red and blue curves are 
almost indistinguishable and correspond 
to the maxima of Cp and Kt respec- 
tively. Near the critical pressure, these 
merge into a single line. 



example the T = 1.8 isotherm plotted in green in fig. ([7]), the two divergences 
of R come closer to each other, and at Tc = 2, they merge at = 1. Beyond 
T = 1, \R\ (~ shows two maxima, as can be seen from the blue curve in 
fig-© which corresponds to an isotherm T = 2.1. The Widom line therefore has 
to be understood in more details here. 

For this, beyond Tc, we study the maxima of \R\ by the equations 



d\R\\ n f9\R\\ n fd\R 



dT Jy ' \dV Jp ' \dT y p 



(33) 



We find that whereas the maxima of \R\ obtained from the first of these relations 
(i.e isochoric maxima of \R\ with respect to the volume) is virtually indistin- 
guishable from the Kt maxima line, those obtained from the second and third 
relations of eq.( l33|l indicate two local maxima for \R\ near criticality. The values 
of V obtained from both these relations (where \R\ maximises) are indistinguish- 
able. Fig. ([9]) summarizes our results. Here, the red curve indicate the locus of 
maxima of Cp, and the blue one is the corresponding quantity for K^. The so- 
lution of the first relation of eq. fl33|) is identical to the blue curve. The solutions 
of the other two relations in eq. fl33|) give the brown and black curves. 

As in the case of the CW model of the previous section, instead of looking 
at isochores, it is more useful to look at the behavior of \R\ with respect to 
a thermodynamic intensive variable, such as the temperature or the pressure. 
In fig. (fTO|) . we have shown the behavior of isothermal \R\ as a function of the 
pressure, for various values of the temperature, with volume being the parameter 
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in the plot. Consider, for example, the red curves in fig.( !TO|) . These are isothermal 
plots of |-R| as a function of the pressure, for T = 1.6, in the two physical domains 
of V (as we have discussed) and correspond to the red curve in fig. ([8]). One of the 
red curves here (the one diverging at P ~ 12.8) corresponds to V < 0.7899 and 
the other red curve is for V > 1.4844. \R\ for the physical branches is thus seen to 
cross at P ~ 12.85. This is thus the value of P for which the correlation lengths 
of the coexisting phases become equal, and we interpret this as the pressure for 
which a first-order liquid-liquid phase transition occurs at T = 1.6. Similar values 
of the pressure can be calculated using isotherms for different values of T. The 
collection of all such points in the T — P plane is the phase co-existence curve 
for the model. 

As we approach criticality, the crossing point of R is pushed to infinity. Be- 
yond criticality, there is no crossing, i.e a single phase exists, but |P| shows two 
maxima with respect to the pressure, i.e there are two Widom lines that orig- 
inate from the critical point. In fig.l lTTl) . we have plotted the Widom lines for 
this system in the P — T plane. Projected on this plane, the maxima of |P| are 
shown in the green and black curves and the other curves, which are the maxima 
of Cp and become almost indistinguishable. All these curves converge at the 
critical point. 

Before we end, we summarize the main results in this section. Here, we 
have shown that information geometry can be effectively used to establish phase 
behavior in liquid systems, in a simple way. The main inputs that went into 
our calculation is the power law behavior of the internal energy of eq.(l2T]) and 
the coefficients appearing in eqs.f l22|) and ( |23|) . We stress here that although 
for illustrative purpose we have chosen a toy model, once these quantities are 
calculated in any liquid system, our method offers an easy way to determine 
phase behavior, and very importantly, the Widom lines, and may be used for 
experimental predictions. We have also seen an important qualitative difference 
between magnetic and liquid systems, namely that in the latter, there might be 
multiple isothermal Widom lines. We note here that this situation is different 
from liquid-gas systems as well, where there is a single locus of maxima of the 
correlation length ^ for isotherms 0. 

5 Conclusions and Discussions 

In this paper, we have studied in details information geometry of magnetic and 
liquid systems. We first established the idea that equality of correlation lengths 
of co-existing phases, calculated via the scalar curvature of the equilibrium ther- 
modynamic state space, is indicative of first order phase behavior in mean-field 
magnetic systems, where the curvature diverges appropriately at criticality. We 
then applied this to liquid systems and predicted liquid-liquid phase transition in 
a simple model. Our main conclusion here is that geometric techniques provide a 
universal new method of characterizing first and second order phase transitions, 
and can also be used to predict the behavior of the system beyond criticality, via 
the Widom line. Our results show that the definition of the latter as the locus of 
maxima of the correlation length is somewhat ambiguous in liquid systems, and 



14 



can lead to multiple lines, which originate from criticality. This is due to the T^/^ 
behavior of the internal energy. This is to be contrasted with liquid-gas systems 
where the prediction of isobaric or isothermal Widom lines are unique (For 
magnetic systems, such multiple locus of maxima occur instead for the specific 
heats). For the 1-D Ising model, we have seen that the Widom line does not 
originate from the critical point. This however, can be attributed to the some- 
what unphysical nature of the model, and this feature is not seen for any of the 
mean-field theories analysed in this paper. It is an important question whether 
the Widom line is unique, and future experiments will probably indicate which of 
the multiple lines is chosen by the system to distinguish between phases, beyond 
criticality. 

In hquid-gas systems, prediction of first order co-existence via the equality of 
the correlation length originates from the idea of Widom \TT\ that near a first or- 
der transition, density fluctuations in one phase of a fluid results in the formation 
of a second phase. The thickness of the interface between the two phases is then 
interpreted as the correlation length, which therefore should be equal measured 
in either phase [6]. Our results indicate that a similar phenomenon happens in 
magnetic and liquid systems as well, and this should be studied further. 

Compared to the more sophisticated computer simulation techniques to study 
liquid-liquid phase transitions advocated over the last two decades, ours is a 
simple method which can be applied to any theoretical equation of state, or 
equivalently, experimental data. Admittedly, we have chosen a simple scenario 
to illustrate our method, but this can be easily generalised to more realistic 
situations, and the qualitative details should remain unchanged. 

It will be of interest to apply our technique to phenomenological models of 
liquid-liquid phase transitions. We note here that the predictions using the scalar 
curvature remain valid as long as the numerical value of \R\ is more than typical 
molecular volumes. For liquid-gas systems, these restrict the use of geometric 
methods beyond a certain range of temperature and pressure, though these are 
typically away from the scaling region. A similar analysis in liquid systems would 
involve analysing data from molecular dynamics simulations. We leave this for a 
future publication. 
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